Homoclinic chaos and energy condition violation 
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In this letter we discuss the connection between so-called homoclinic chaos and the violation of 
energy conditions in locally rotationally symmetric Bianchi type IX models, where the matter is 
assumed to be non-tilted dust and a positive cosmological constant. We show that homoclinic chaos 
in these models is an artifact of unphysical assumptions: it requires that there exist solutions with 
positive matter energy density p > that evolve through the singularity and beyond as solutions 
with negative matter energy density p < 0. Homoclinic chaos is absent when it is assumed that 
the dust particles always retain their positive mass. In addition, we discuss more general models: 
for solutions that are not locally rotionally symmetric we demonstrate that the construction of 
extensions through the singularity, which is required for homoclinic chaos, is not possible in general. 

PACS numbers: 04.20.-q, 98.80.Jk, 98.80.Cq, 04.20.Ha 



The dynamics of locally rotationally symmetric (LRS) 
and diagonal Bianchi type IX models with a perfect fluid, 
usually chosen to be dust, and a positive cosmological 
constant has been discussed recently in a series of arti- 
cles 0, U H • It was claimed that Einstein's static uni- 
verse is a chaotic scatterer so that "homoclinic chaos" 
is generated in the models. In another article [f|, LRS 
type IX models with general equations of state (compris- 
ing dust and a positive cosmological constant as a special 
case) were re-investigated by employing a scale-invariant 
formulation. This permitted the use of rigorous dynam- 
ical systems methods which led to results that excluded 
the possibility of homoclinic chaos; the dynamics of the 
models was found to be predictable. In spite of these 
rigorous results, claims about homoclinic chaos persist 
in the literature. The purpose of this letter is to resolve 
these contradictions. 

In the first part of this letter we show that homoclinic 
chaos in LRS Bianchi type IX models is an artifact of 
unphysical assumptions: for the models to exhibit homo- 
clinic chaos solutions have to be extended through the 
singularity where it must be assumed that the mass of the 
dust particles changes sign. The assumption of standard 
energy conditions (such as p > for dust) excludes homo- 
clinic chaos and leads to the results of |5|. In the second 
part, in addition to the issues of energy condition viola- 
tion, we point out that the feasibility of the (unphysical) 
construction leading to homoclinic chaos entirely relies 
on the assumption of LRS symmetry, i.e., extensions of 
solutions through the singularity are not possible in more 
general models — the special structure of LRS type IX 
models is misleading. 

For LRS Bianchi type IX models the metric can be 
written as 



ds 2 



-dt 2 + A 2 (t) {uj 1 ) 2 + B 2 {t) [(a, 2 ) 2 + (w 3 ) 2 ] , (1) 



where the uj 1 are invariant Bianchi type IX one-forms. 
We consider models whose matter content is non-tilted 



dust with energy density p, and a positive cosmological 
constant A. We set c = 1 = 8irG. The Hamiltonian for 
this class of models reads 
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where pa, Pb are the momenta conjugate to A, B. 
The constant Eq corresponds to a matter term, Eq = 



2pAB 2 = const, 
Bianchi identities, 
by 
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which follows from the contracted 
Einstein's field equations are given 
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The system (0 admits Einstein's static universe as a 
solution, represented by the fixed point E: pa = Pb = 0, 
A = B = I/VAX, E = 1/VaK- In 0, the field equations 
were reformulated in terms of scale-invariant variables in 
order to obtain a regular system. This allowed the use of 
rigorous dynamical systems methods and thus resulted 
in a comprehensive description of the global dynamics; 
in particular, the dynamics of models initially close to E 
was found to be predictable and non-chaotic. Since the 
results are completely invariant and do not rely on the 
particular choice of variables, the findings oflal refute 
earlier claims about "homoclinic chaos" [l|, |3, y, |j] . 

The claims about homoclinic chaos rely on the asser- 
tion that there exist infinitely many solutions of J2J that 
evolve from an initial state close to E, leave any neigh- 
borhood of E, and return to a final state close to E again. 
It is argued that this behavior is the basis of chaos in the 
dynamics because of the chaotic distribution of "homo- 
clinic" initial data points among initial data close to E. 
In the following we demonstrate that the asserted "homo- 
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clinic solutions" can only exist under unphysical assump- 
tions. For pedagogical reasons, in our considerations we 
focus on the most recent publication of 0,11110, i.e., ||; 
in particular we adopt the notation and terminology. 1 

The analysis of 4] is crucially based on the observation 
that the system @ is regular for all A including A = 0. 
This prompts the authors to not restrict the domain of 
the variable A to A > ("truncation at A = is bound 
to a loss of information" ) , but to consider the regions 
A < and A = as a part of phase space ("for the 
Hamiltonian dynamics [the region A < 0] constitutes an 
essential part of the phase space"). 

As a first comment, let us note that the regularity of 
the equations J2J at A = is inseparably connected with 
the assumption of LRS symmetry. Generalizing JJJ to the 
full diagonal type IX case with spatial metric A 2 ^ 1 ) 2 + 
B 2 (cu 2 ) 2 + C 2 (lu 3 ) 2 yields 
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for the kinetic term of the Hamiltonian. (For B = C, 
Pb = Pc, the kinetic term of J5J is recovered by noting 
that the new momentum ps differs from the original one 
by a factor of 2.) We observe that the non-LRS Hamil- 
tonian and the resulting equations are singular at A = 
(and B = 0, C = 0). Any deviation from LRS thus in- 
validates the statement that "the regions of phase space 
A > and A < join smoothly in A = 0." This will be 
further discussed below. 

In general relativity, for LRS line elements of the 
form one usually makes the restriction A > (and 
B > 0); in particular, Hamiltonian approaches typically 
use exponential representations of the scale factors, e.g., 
A = e , see, e.g., EH US OH. it is tacitly, and cor- 
rectly, assumed that the restriction to positive scale fac- 
tors is no loss of generality. This is due to the existence 
of discrete symmetries: the system © is invariant under 
(A,p A ,t) h-> ~{A,p A ,t) and {B,p B ) t-> ~{B,p B ) (where 
the latter reflects the invariance of the non-LRS sys- 
tem under the simultaneous transformations (B,pB,t) i— > 
— (B,pB,t) and (C,pc,t) t— > — (C,pc,t)). As a conse- 
quence, solutions with A < arise from solutions with 
A > via this symmetry map; in this context it is im- 
portant to note that Eq h- > —Eq and that p is invariant 
under the transformations. 

The phase space considered in |4| is the union of the 
regions A > and A < 0, "smoothly joined on A = 0" ; 
in particular, A = is not regarded as a singularity from 
a Hamiltonian point of view. The authors give numerical 
evidence for the existence of an infinite set of solutions 



1 We also adopt the term "homoclinic solutions" here, although 
this is a deviation from standard dynamical systems nomencla- 
ture. 



with initial data close to the fixed point E (i.e., with 
A > 0) that pass through ^4 = 0, enter the region A < 0, 
re-enter A > through A = 0, and eventually enter a 
neighborhood of E again. It is the existence of these 
"homoclinic solutions" that is argued to cause "homo- 
clinic chaos" in the model. For the homoclinic solutions 
Eq = const > holds; in fact, Eq is close to l/v4A, 
which is the value characterizing Einstein's universe E. 

The positivity of Eq has far reaching consequences: 
from Eq = 2AB 2 p we obtain 
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Eq 



2AB 2 ' 



(4) 



which implies that p < 0, if A < (since Eq > 0). 
Therefore models described by solutions with A < vio- 
late the (weak/null/strong/dominant) energy condition; 
indeed, since the matter content is dust, p < implies 
that the dust particles have negative inertial and gravi- 
tational mass. The homoclinic solutions described in 
evolve between the regions A > and A < 0, hence p > 
changes to p < (and back) during the evolution. How- 
ever, the change of the density of the dust particles from 
positive to negative is not continuous: when A \ we 
find p — > oo, when A /" we have p — > — oo. (This is 
not immediate from (J2J|, but follows from the fact that 
the spatial volume density AB 2 converges to zero as the 
singularity is approached.) Hence extending a solution 
from A > via A = to A < (and back) corresponds 
to joining p = +oo and p = — oo (and again — oo with 
+oo). 

In jj] the authors emphasize that it is essential to in- 
corporate the regions A < and A = into the analysis, 
since the "dynamical phenomena that occur in A < 
are present in A > transported there via the Hamil- 
tonian flow" . Naturally, there exist numerous examples 
in physics where the consideration of unphysical states is 
used as an abstract mathematical tool that makes possi- 
ble or facilitates the analysis of physical states. Here the 
case is totally different. The system © is an autonomous 
system, which makes the flow a local phenomenon. In 
particular the phase space picture in the region A > 
depends exclusively on A > 0; it is a basic principle that 
(changes of) the flow in A < cannot affect the phase 
space picture in A > 0. 

Though irrelevant, the assertion that the system © 
in the region A < is associated with complex dynamics 
(see, e.g., Figs. 17 and 18 in Q) is likely to be correct, but 
considering that we are dealing with dust particles with 
negative masses this should perhaps not be too surpris- 
ing. Note that the occurrence of negative masses is not 
primarily due to A < 0, but due to the fact that A < 
while Eq is kept positive. Recall in this context that 
application of the symmetry map (A, Eq) h- > — (A,Eq), 
cf. the discussion above, would leave the dynamics (and 
the positivity of p) unchanged. 

We conclude that there exist two attitudes one can 



3 



take to home-clinic chaos in LRS type IX models with 
dust and a positive cosmological constant. 

(i) A conservative attitude, which is the prevailing 
opinion of the classical general relativity commu- 
nity; this is also the view we hold. Restricting the 
phase space to positive scale factors is no loss of 
generality; A = represents a singularity of the 
model, since p — > oo when A — ► 0. Extending 
a cosmological model as a solution with A < 
(and Eq > 0) amounts to joining p = +oo with 
p = — oo, and it must be assumed that the dust 
particles replace their originally positive mass by 
negative mass. Any such extension is unreasonable 
from the point of view of physics, and irrelevant for 
the description of physical states from the point of 
view of mathematics. 

(ii) An unconventional attitude. The phase space is 
given by the union of A > 0, A = 0, and A < 0; A = 
is not a singularity in the sense that the evolution 
of cosmological models does not stop. There exists 
a mechanism that transports the model (which is 
characterized by Eq > and thus has positive mass 
originally) through A = (p — ±oo) to a model 
with negative mass, and back. The constructed 
spacetimes are of the form R x M, where the man- 
ifold M is Riemannian except for at a discrete set 
of times, where the metric on M is degenerate. At 
times of degeneracy, the density of the dust par- 
ticles changes sign from p > via p = +oo and 
p = — oo to p < (or reversed). 

If one takes view (i), it can be proved, see that 
no homoclinic chaos exists and that the dynamics is pre- 
dictable. The results in 5] are completely invariant and 
do not rely on the particular scale-invariant formulation; 
however, the results do rely on the assumption that the 
dust particles retain their positive mass. 2 

View (ii) is associated with severe problems. First, we 
note that although formally, i.e., from a purely mathe- 
matical point of view, the extension of the scale factor A 
through A = is feasible, this extension relies crucially 
on the assumption of LRS symmetry. When LRS sym- 
metry is broken (as it certainly is for generic models) it is 
far from clear whether an extension of generic solutions 
through the singularity exists. Let us elaborate on this 
issue. 



2 The essence of view (i) pervades all theories that model phenom- 
ena with differential equations: although the equations might be 
well-defined (and regular) for a range of a variable that is greater 
than its physical range, the physicality condition prohibits the 
extension of solutions. Examples abound not only in physics; for 
instance, in population models in mathematical biology llll . so- 
lutions will never be extended to negative numbers of individuals 
of a species. 



In the LRS case, the regularity of (I2J and |3J in A al- 
lowed |J| to also regularize the problem in B by a change 
of the time variable that only involved B. In this for- 
mulation, B = was an invariant subset, so that orbits 
could not pass through B = 0; in contrast, the set A = 
was not invariant, which was exploited to extend orbits 
beyond the singularity A — 0. In the full diagonal type 
IX case the equations are singular at A = 0, B = 0, and 
C = 0. It is possible to regularize the equations by a 
change of the time variable (in a spirit similar to the B 
regularization in 4]). However, this re-parametrization 
must necessarily involve all three scale factor due to per- 
mutation symmetry, e.g., one could choose a lapse pro- 
portional to ABC. As a consequence, A = 0, B = 0, 
and C = become invariant subsets, which prevents so- 
lutions from changing the sign of the scale factors (and 
thus from changing the sign of p). Hence, in the non-LRS 
case, an extension of solutions through the singularity 
lacks a mathematical foundation. 

That LRS models are of a special nature is well known. 
In particular, singularities that occur in the LRS case 
are very different from those occurring in the diagonal 
case. In the LRS case solutions generically approach the 
singularity along a so-called Taub asymptote. This sug- 
gests that the singularity may be a weak null singularity, 
see ^3 an< i references therein, and that is possible to ex- 
tend the spacetime through the singularity in a C° man- 
ner to the Minkowski spacetime. It is thus not inconceiv- 
able that the metric can be extended, but it is unclear 
whether such a construction bears any relation to the 
extension of the scale factor A critically discussed above. 
The problem of whether the extension of A leads to an 
extension of the metric is not discussed in ; clearly, the 
extension of A is completely irrelevant unless it leads to 
an extension of the metric. 

In the diagonal (non-LRS) case, on the other hand, 
the situation is completely different: solutions generi- 
cally approach the singularity in a Mixmaster fashion 
and the singularity is a spacelike scalar curvature singu- 
larity 13]; no extension of the metric is possible. Since 
the Mixmaster behavior is typical for generic cosmologi- 
cal singularities, see e.g. 0] and references therein, these 
considerations seem to rule out the type of manipulations 
and arguments presented in |4j for generic models. It ap- 
pears that the special LRS case is quite misleading when 
it comes to spacetimes and phase space extensions. 

Even in the LRS case, not all solutions possess ex- 
tensions of the kind discussed above; a Taub asymp- 
tote is a prerequisite. This is because for such solutions 
B y4 when the singularity is approached; solutions with 
B — ► 0, e.g. solutions with an isotropic singularity, are 
inextendible. This creates a link to the related question 
of why it is possible for unphysical solutions (p < 0) to ex- 
pand from a singularity, and then rccollapse to a singular- 
ity at A = (and subsequently re-enter the physical part 
of the phase space) . At first this is quite surprising since 
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dust particles with negative mass generate anti-gravity, 
and one would thus expect such solutions to expand to 
a state infinite dilution, p /* 0, which is confirmed by 
numerics for typical solutions. However, the solutions 
found numerically in have Taub asymptotes and are 
thus vacuum dominated with extreme shear; since shear 
represents gravity generating gravity, these solutions can 
recollapse to the singularity and enter the physical region 
again. Note, however, that this scenario is only correct 
for solutions with Taub asymptotes with initial data such 
that the shear dominates during a sufficiently long time. 

Finally, let us point out that, irrespective of the math- 
ematical considerations, at present there does not exist 
any physical theory that could explain the traversing of 
the singularity in the way implicitly suggested by adher- 
ents of view (ii). Note that certain speculative mech- 
anisms like "quantum tunneling of the wave function of 
the universe" (see 0|) are insufficient, since they lack an 
explanation for the metamorphosis of the dust particles 
from positive to negative mass. In addition, in view of 
the discussion above, such a mechanism is a serious candi- 
date only if it can explain the behavior of generic models 
and not merely the exceptional LRS type IX models. 

It seems to us that the burden of making a convincing 
case for homoclinic chaos in Bianchi type IX models lies 
with adherents to view (ii); considering the above, this is 
a formidable challenge indeed. 

Note added. 

After submission of this letter an erratum to was 
published, see [J^. In this small note the authors state 
that the physical interpretation on the existence of ho- 
moclinic chaos has to be corrected, but the authors do 
not give any explanation for this statement (see, how- 
ever, the first part of the present letter); furthermore, it 
is not discussed that the corrected interpretation also af- 
fects the status of homoclinic chaos in non-LRS models 



(see, however, the second part of the present letter). 
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